Quantum Monte Carlo study of the visibility of one-dimensional Bose-Fermi mixtures 
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The study of ultracold optically trapped atoms has opened new vistas in the physics of correlated 
quantum systems. Much attention has now turned to mixtures of bosonic and fermionic atoms. 
A central puzzle is the disagreement between the experimental observation of a reduced bosonic 
visibility Vb, and quantum Monte Carlo (QMC) calculations which show Vb increasing. In this 
paper, we present QMC simulations which evaluate the density profiles and V;, of mixtures of bosons 
and fermions in one-dimensional optical lattices. We resolve the discrepancy between theory and 
experiment by identifying parameter regimes where Vb is reduced, and where it is increased. We 
present a simple qualitative picture of the different response to the fermion admixture in terms of 
the superfluid and Mott-insulating domains before and after the fermions are included. Finally, we 
show that Vi, exhibits kinks which are tied to the domain evolution present in the pure case, and 
also additional structure arising from the formation of boson-fermion molecules, a prediction for 
future experiments. 

PACS numbers: 03.75.Mn, 03.65.Yz, 03.75.Hh, 71.10.Pm 



It has been widely suggested that the strong correla- 
tions responsible for magnetism, superconductivity, and 
the metal-insulator transition in the solid state can be 
studied via ultracold optically trapped atoms. Indeed, 
this idea has been successfully realized in the context of 
both bosonic and fermionic atoms. In the former case, 
the transition between condensed (superfluid) and insu- 
lating phases was demonstrated through the evolution 
of the interference pattern after the release and expan- 
sion of the gas JM. Initial studies focused on the height 
[l[ and width 0] of the central interference peak, with 
later work looking at the visibility V, which measures 
the difference between the maxima and minima of the 
momentum distribution function S(k) 0, 0,01 ■ Interest- 
ing "kinks" are observed in V which are associated with 
the redistribution of the density as the superfluid shells 
evolve into insulating regions 0, ■ For trapped fermions 
BiMinLGl, Mott phases could also form [3 El 13 > 
however, without a signal in 5(k). Instead, the evolution 
of the kinetic energy has been proposed as a means to 
pinpoint the transition [3- 

Attention has turned at present to multicomponent 
systems, which offer a rich set of phenomena includ- 
ing Bose-Einstein condensation (BEC)-BCS crossover for 
two attractive fermionic species and Fulde-Ferrell-Larkin- 
Ovchinnikov phases in situations with two imbalanced 
fermion populations. Two recent experimental papers re- 
port the effect of the addition of fermionic 40 K atoms on 
the visibility of bosonic 87 Rb in a three-dimensional trap 
[13, EH- The basic result is a decrease in the bosonic 
visibility Vb driven by the fermion admixture. A large 
number of qualitative explanations has been put forth 
for this phenomenon, including the localization of the 
bosons by the random fermionic impurities, the segmen- 
tation of the bosonic superfluid, the adiabatic heating of 
the bosonic cloud when the lattice depth is increased in 
the presence of the two species, an enhanced bosonic mass 



due to the coupling to the fermions, and the growth of 
Mott-insulating regions. A fundamental difficulty is that 
exact quantum Monte Carlo (QMC) calculations show 
an increase in Vb [l9j . in disagreement with experiment. 
In this paper we resolve this issue. 

The behavior of Bose-Fermi mixtures has attracted 
considerable theoretical attention. The Hamiltonian was 
derived and its parameters linked to experimental quanti- 
ties by Albus et al. [2~oj . The equilibrium phase diagram 
has been studied using mean-field theory and Gutzwillcr 
decoupling [2(| HH, [HJ , perturbation theory [2l|, dynami- 
cal mean- field theory (DMFT) [3 , exact diagonalization 
[23| . quantum Monte Carlo methods [H, IS, ll , and 
density matrix renormalization group (DMRG) [3 [27j • 
The results of these studies include the observation of 
Mott-insulating phases at "double half-filling" , anticor- 
related winding of the two species of quantum particles, 
molecule formation, and precise determination of the ex- 
ponents characterizing correlation function decay in the 
different phases. The behavior of the visibility was ad- 
dressed by Pollet et al. [3, who find interesting non- 
monotonic structures with fermion density. However, Vb 
is always increased relative to the pure case f28j . 

In this paper we explore the visibility of Bose-Fermi 
mixtures in one dimension using QMC simulations with 
the canonical worm algorithm gMHl- While previ- 
ous QMC studies have reported a growth of Vb, we show 
that a significant reduction, such as seen experimentally, 
is also possible without invoking temperature effects [l9( . 
The enhancement (reduction) of Vb caused by the disrup- 
tion (inducement) of the bosonic Mott-insulator phase 
by the boson-fermion interactions. Vb also exhibits kinks 
reminiscent of those in the pure boson case. In the sub- 
sequent sections we write down the Hamiltonian and ob- 
servables and briefly discuss the QMC algorithm. We 
then present the evolution of Vb with fermion concentra- 
tion, its interpretation in terms of the bosonic density 
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FIG. 1: (Color online) Comparison of the density profiles at 
Ubb = S.3t and Ubf = —5.0t for Nf = 0,3,5 fermions on an 
80-site chain with 40 bosons. The Mott insulator at the trap 
center for the pure bosonic case is destroyed by the addition 
of fermions. This drives the increase in the visibility. 



profiles, and evidence for the formation of a molecular 
superfluid in the trap center. 

The one-dimensional (ID) Bosc-Fermi Hubbard 
Hamiltonian is 1201 



H 
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where bj (Sj) and dj (cj) are the annihilation (creation) 
operators of the bosons and (spin-polarized) fermions at 
lattice site j, respectively, and n 



'/ 



= cjc, 



are the corresponding number operators. The first two 
terms of Eq. ((T|) describe bosonic and fcrmionic nearest- 
neighbor hopping. The curvature of the trap is W , and 
the coordinate of the jth site is given by Xj = ja, where a 
is the lattice constant. Ubb and Ubf are the on-site boson- 
boson and boson-fermion interactions. In this work we 
consider 80-sitc chains with the nearest-neighbor hopping 
set to be identical for fermions and bosons (f & = tt = t = 
1) and trapping potential W = O.Olt. 

In the canonical worm algorithm [29|, H3, HH employed 
in our calculation, operator expectation values are sam- 
pled through the introduction of open-ended world lines 
that extend over equal imaginary time into a path inte- 
gral expression for the partition function. The properties 
we study include the kinetic, potential, and trap energies, 
the density profiles, and the visibility, 



V 



5 a 



S n 



(2) 



where SVax and S m i n are the maximum and minimum 
values of momentum distribution function, 



s(k) = i^y fcfe - 

3,1 
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FIG. 2: (Color online) (a) Bosonic visibility Vb as a function 
of the number of bosons for Nf = and 9 fermions with fixed 
Ubb = 12. Of and Ubf — —5.0t. (b) Bosonic density profiles 
for Nb — 26 bosons and Nf = 0, 9 fermions. The addition of 
the fermions induces Mott-insulating behavior in the bosons. 
The key consequence is a decrease in V;, for Nf — 9 relative 
to Nf = 0, similar to that seen in the experiments. 



The enhanced visibilities with fermion concentration 
reported previously [l9j are in contrast with the trend to 
reduced Vb measured experimentally [UGH- ^ n Fig- HI 
we see the origin of this effect in a system with 40 bosons: 
the visibility enhancement at large Ubb is caused by the 
destruction of the Mott phase at the trap center by the 
fermions. It is natural to conjecture that if < 1 
at the trap center the additional attraction due to the 
fermions could induce Mott-insulating behavior and re- 
duce V. In Fig. HJa), we show that this expectation is 
correct. Here, we fix Ubb = 12i and increase Nb for both 
the pure case and for a system with fermion number fixed 
at Nf = 9 and boson-fermion interaction at Ubf = —5t. 
What we observe is that in a window where the boson 
central density is approaching n\ = 1 the bosonic visi- 
bility is decreased by the presence of the fermions. The 
cause is clear: if the bosons arc poised just below Mott- 
insulating behavior, then the fermions can induce it. This 
is supported by a comparison of the density profiles in 
Fig. [lb). 

The primary mechanism through which fermions affect 
Vb is the local adjustment of the site energy and hence 
of the local bosonic density. This is an effect which oc- 
curs regardless of the dimensionality. Hence, we expect 
aspects of our conclusions to be relevant to experiments 
in higher dimension [l^|. While we have shown a de- 
creased visibility similar to that seen experimentally, the 
enhancement of visibility may be the more generic be- 
havior in one dimension. In the one-dimensional "state 
diagram" of the purely bosonic case [32| the area of pa- 
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FIG. 3: (Color online) (a) Bosonic and fermionic visibilities 
and bosonic 5 ma x as functions of Ubb/t for a system with 40 
bosons, 3 fermions, and ZTj,/ = —5.0t. The "plateau" regions 
where the rate of reduction of V is reduced are due to freez- 
ing of the density profiles (see text). The fast decrease after 
Ubb/t ~ 9.3 is due to the formation of a Mott-insulating re- 
gion in the central core, which is fully formed and indicated 
by the arrow at Ubt/t — 9.6. (b) Boson density profiles at five 
different values of Ubb/t. 



rameter space occupied by the phase with a Mott plateau 
of n b = 2 is very narrow. Thus the prospect for the 
fermions to drive the system into this phase is limited. 

In the case of a pure bosonic system 0] , the change in 
visibility with the boson-boson interaction strength Ubb 
is not smooth, but is accompanied by "kinks." These 
kinks are associated with a freezing of the density profile 
when the transfer of the bosonic density from the trap 
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FIG. 4: (Color online) (a) Trapping (Strap) and interaction 
(-Eint) energies as functions of Ubb/t for the system of Fig. [3] 
(b) Bosonic (I?£ in ) and fermionic (-Ejfin) kinetic energies. 



40 - 




9 10 11 

U b b/t 

FIG. 5: (Color online) Comparison of (a) bosonic visibility 
Vb, (b) the trapping energy, and (c) the interaction energy for 
Nf = 0, 3, and 5 fermions on an 80-site chain with 40 bosons 
and fixed Ubf = — 5.0t. The arrows in panels (a), (c), and 
(b), respectively, denote the locations of the kinks (the onset 
of rapid change in energy and visibility) for Nf = 0, 3, and 5 
fermions. 



center is interrupted by the formation of Mott insulator 
shoulders. In Fig. [3j the behavior of the visibilities and 
density profiles with Ubb in the presence of fermions is 
shown. Vb decreases with Ubb as the interactions reduce 
the quasicondensate fraction S^^. Conversely, the inter- 
actions enhance S^ax and V/ increases with Ubb- There 
are, however, additional kinks in the case when fermions 
are present whose origin we shall discuss below. Figured] 
helps to quantify this freezing by showing the evolution 
of the trap, interaction, and kinetic energies with Ubb- 
These energies exhibit a sequence of plateaus and rapid 
drops corresponding to the kink locations in Fig. [31 

Figure [5ja) compares the visibility evolution for the 
pure bosonic case (Nf = 0) with two different fermion 
numbers Nf = 3 and 5. For Nf = 3, the kink at lowest 
Ubb = 6. It is associated with the initial emergence of the 
Mott shoulders. This kink coincides with one in the pure 
bosonic case Nf = because the shoulders form outside 
the regions occupied by the fermions at the trap center. 
For Nf = 5, the width of the fermion density is compa- 
rable to the size of the bosonic supcrfhiid in the center 
of the trap, and the kink at Ubb = 6. Of signifies a freez- 
ing of the bosonic density but not the formation of the 
Mott shoulders. Instead, the kink visible at Ubb ~ 6.8t 
is responsible for the formation of the Mott shoulders. 
This shift to higher Ubb is expected since the attractive 
Ub f delays the transfer of bosonic density out of the cen- 
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FIG. 6: (Color online) Bosonic and fermionic density profiles 
for (a) Ubb = 6.3t, (b) U bb = 7.2t, (c) U bb = 8M, and (d) 
Ubb = 9.5t with 5 fermions, 40 bosons, and fixed U b f = —5. Of. 
The fermionic density is offset and the dashed gray line in- 
dicates zero density. The densities match in the center of 
the trap, with the region of coincidence decreasing as U bb in- 
creases. 



tcr. Figures [5^b) and EJc) compare the components of 
the energy. Each plateau signifies that the bosonic and 
fermionic densities are frozen over the range in Ubb- Past 
experiments fl7l . [l8| did not have the resolution to ex- 
hibit the kinks we have seen in these simulations; how- 
ever, they might be seen with improved accuracy. 

We also note in Fig. \5\ that the number of plateaus is 
directly related to the number of fermions in the Bosc- 
Fermi mixture and that each plateau is roughly the same 
size, indicating that bound pairs of bosons and fermions 



are being destroyed as Ubb is increased. This conclusion 
is substantiated in Fig. [TjJ where we compare the bosonic 
density with the fermionic density. The fermionic density 
is offset by a constant to emphasize the near perfect over- 
lap in the densities near the center of the trap, indicating 
that the trap center is populated by a molecular super- 
fluid (MSF). fndeed, at the weakest coupling, Fig. [HJa), 
the fcrmion density precisely equals the excess boson den- 
sity above the commensurate Mott value ni^ = 1. When 
Ubb is increased, moving from one plateau to another in 
Fig. [5j the MSF region in the center of the trap shrinks. 
For the kink at highest U b b (~ 9.8i for N f = 5), the 
MSF region is destroyed, the bosonic density is a Mott 
insulator, and the fermionic visibility V/ — » 1. 

In summary, we have shown that the visibility of Bose- 
Fermi mixtures can be enhanced or reduced by the boson- 
fermion interactions depending on whether the bosonic 
density in the pure case is above or below commensu- 
ration. This result resolves a fundamental disagreement 
between experiment and QMC simulations. There arc 
numerous kinks in the visibility and the different en- 
ergies that result from freezing of the density profiles. 
While our bosonic component is sufficiently large so that 
our results are converged with respect to lattice size, the 
number of fermions is much smaller. It is possible that 
the kinks will merge together and be less easy to observe 
in a larger system. The density profiles near the kinks 
show direct evidence for a molecular supcrfluid in the 
center of the trap and that a larger Ubb is required to de- 
stroy the bound pairs with larger Nf and subsequently 
induce Mott-insulating behavior. 
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